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Abstract. Transient phenomena in quantum mechanics have been of interest to one of the 
authors (MM) since long ago and, in this paper, we focus on the problem of a potential V- 
which for negative times gives rise to bound states and is suddenly changed at t = to 
a potential V+ which includes V- plus a perturbed term. An example will be the deuteron 
(where the proton and neutron are assumed to interact through an oscillator potential) 
submitted to a sudden electrostatic field. The analysis for t > can be carried out with the 
^ ' help of appropriate Feynmann propagators and we arrive at the result that the separation 

CsJ , between the nucleons has an amplitude that depends on the intensity of the electrostatic 

field, but its period continues to be related with the inverse of the frequency of the oscillator 
proposed for the interaction. A general approximate procedure for arbitrary problems of 
this type is also presented at the end. 
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^3 . 1 Introduction 

The word "Transient" will mean for us here that the expectation values of measurable observable 
^ , or transition probabilities in quantum mechanics will be time dependent, in contrast with the 

^ I "stationary" case in which they are not. Prom the very origin of quantum mechanics, in the 

ninety twenties, these transient phenomena could be discussed through the time dependent 
versions of the equations of motion in the Heisenberg or Schrodinger picture. 

The most natural way of discussing the transient phenomena is through the version of quan- 
tum mechanics that Feynman ^ proposed in the ninety forties. For him the fundamental 
concept was a propagator, also sometimes called Feynmans kernel, that relates an initial state 
characterized by the wave function ip{x' , t') (where x' is a notation for all the coordinates of the 
system and t' its time) with the final state 'ijj{x,t) with t > t' . This propagator was denoted by 
K{x,t;x',t') and Feynman gives a procedure to calculate it starting from a Lagrangian formu- 
lation (in terms of generalized coordinates and velocities) of the classical problem and summing 
all paths suggested by the action principle between x' and x. If the classical Lagrangian, is not 
explicitly dependent on time, the propagator only depends on the difference of the times i.e. 
t — t' so, without loss of generality, we could take t' = and write 



K{x, t; x', 0) = K{x, t; x'), i/j{x', 0) = Tp{x') 
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so that we have 



tp{x,t) = K{x,t,x')'tp{x')dx' 



(1) 



which imphes that the propagator K{x,t,x') satisfies a time dependent Schrodinger equation 
with the initial condition 



The use of the propagator concept has been enhanced by the appearance a few years ago 
of a "Handbook of Feynman Path Integrals" [2, and the power of the concept, in the case of 
a one dimensional free particle, is illustrated in the Appendix for a compact derivation of the 
problem of "Diffraction in time" , which one of the authors (MM) |H] analyzed long ago by other 
methods. 

The general type of problems that we wish to discuss here is that of the potential V-{x) 
which for negative times t < is one that admits bound states and that at i = is suddenly 
changed to a potential V+{x) that contains V-{x) plus some type of external interaction given 
by V+{x) — y_ (x). If our vector x is given in cartesian coordinates the kinetic energy is a sum of 
terms (m^x^g) where q = 1, . . . , m is the index for the dimension m of our single particle vector 
space and s = 1, . . . ,n the index for the number of particles. Thus our classical Lagrangian is 
well defined and we can look through the tables of reference [2] to see whether the propagator 
is available. We thus proceed to analyze the following problem. 

2 The deuteron subjected to a sudden electrostatic field 

The first physical problem that came to our attention along the lines of the last paragraph of the 
previous section, was the hydrogen atom in a sudden constant electrostatic field. Unfortunately 
the propagator for the Coulomb problem alone is already quite complicated, and more so if 
a linear electrostatic potential is added. 

We thus wanted to focus on problems in which the propagator can be fully and simply 
determined, and that led us to Lagrangians that are quadratic in the coordinate and velocities 
variables, of which there is an extensive list of propagators in reference 0. 

Why are quadratic Lagrangian (which imply also quadratic Hamiltonians in the coordinate 
and momentum variables) give propagators of the Gaussian type 0? One answer is given 
in a paper by one of the authors (MM) and C. Quesne jlj. The time evolution associated 
with a classical Lagrangian is given by a canonical transformation which conserves a symplectic 
metric. For quadratic Lagrangians or Hamiltonian this canonical transformation is linear and 
it provides the dynamical group of the problem. To translate it to quantum mechanics we have 
to take its unitary representation and this gives a Gaussian propagator [3| . An example of this 
problem is the Lagrangian or Hamiltonian of the harmonic oscillator. 

Thus we were led to the deuteron which we can consider, as is usually done in nuclear physics, 
as a system of neutron and proton, of essentially the same mass m and interacting through an 
harmonic oscillator potential of frequency u>. 

With units in which h = m = c = 1 (c being the velocity of light so that everything will be 
dimensionless) and if we apply suddenly at i = an electrostatic field, in the direction z, acting 
solely on the proton, the Lagrangian will be 



where the indexes 1 and 2 will correspond respectively to the proton and neutron and £ is in 
our units the intensity of a linear potential between the plates of a condenser multiplied by the 
charge of the proton. 



K{x,0;x') = 5{x - x'). 
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We proceed now to consider an orthogonal transformation leading essentially to relative and 
center of mass coordinates 



1 



(ri - r2), 



V2' 

which transforms C to 



1 



1 



R 



8 

T2 



1 

7^ 



+ 



2 



72 



Z 



(2) 



where the components of the vectors r and R are indicated respectively as 

r = ix + jy + kz, R = iX + jY + kZ 

with i, _7, k being unit vectors in the directions indicated. 

If a Lagrangian can be expressed as a sum of two terms depending on different observables 
the propagator is given by the product of the propagators associated with these two terms. 
Furthermore, as the squares of the vectors r, R are given by 



2 2,2,2 
r = X +y + z , 



it is convenient to express our problem in cartesian coordinates where the Lagrangian be- 
comes 



C 



1 



±2 



2 2\ , ^ / -2 

■uj X ) + i^{y 



2 2 

■ UJ Z 



V2£z) 



:X^ + ^Y^ + UZ^-V2£Z) 



(3) 



Thus, from the observation at the beginning of the previous paragraph, the propagator as- 
sociated with C will be the product of the one dimensional propagators associated with the six 
terms appearing in Q. For 



■uj'^x'^) the propagator is given in [21 p. 178, formula (6.2.33)] 



as 



UJ 


^ J 




.27rz sin Lot- 


exp ^ 


~2i 



(x^ + x''^) cot LVt 



2xx' 
sin Lvt 



(4) 



For 

,2^2 



Lo'^y'^) the propagator is again of the form but x, x' replaced by y, y' . For 



\{z^ — uP'z^ — \f2Ez) we complete the square introducing the variable z by 

z = z^ {£/V2io'^) 
and thus the one dimensional Lagrangian becomes 
£^ 



u?z^\ + 



4a;2 



(5) 



(6) 



as z = i. Furthermore for the constant term (£?' jAuP'^ we can apply the relation between 
propagators and Green functions of the time dependent Schrodinger equation, mentioned in p. 2 
of reference [2], to see that it contributes the phase term exp[— ^/4ti;^)i] while the remainder 
in Q is just the one-dimensional Lagrangian of the oscillator for the propagator of which we 
can use equation ((3} replacing x, by z! . It is convenient to express the propagator involving 
the relative coordinates in terms of only barred variables if we add the definitions 



X = X, y = y 
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as in these coordinates we have only the oscillator potential for x and y. With the definitions 

= + + z^, f = ix + jy + kz 

and the previous discussion concerning the phase factor associated with (i5^/4a;^) we obtain then 
that the part of the propagator related to the relative coordinate r becomes 

2r • f'' 



exp [-i{£^/4io^)t\ 



Ijj 


3/2 1 




-27ri sinwi. 


exp < 


~2i 



+ r' ) cot u!t 



sin Lot 



(7) 



Turning now to the expression (j^J with capital letters we start with which from refe^ 



rence |H p. 174, formula (6.2.10)] becomes 



{2mt) 



-exp 



2t 



{X -X 



l\2 



(8) 



For ^Yia same formula ^ but with X, X' replaced hjY ,Y' . For (^Z^ - ^Z) we 



use reference (21 p. 175, formula (6.2.18)] to get 



(27rzt) 



- exp < I 



{Z - Z ') 
2t 



l\2 



JiLiz + z')- — 

2^/2^ ' 48 



(9) 



Multiplying ©, the corresponding expression for y, Y' , and (jUJ we get for the center of mass 
part of the Lagrangian the expression 

'{R-R!f 



(2^it)3/2 



exp < I 



2t 



£t SH^ 
{Z + Z'] 



2y/2 



48 



(10) 



Now the full propagator associated with the Lagrangian Q will be the product of Q and ()1U() . 
For a deuteron at t > the ground state is 

^exp (-^w^^) exp{iK ■ R), (11) 
where A is the normalization of the Gaussian term given by 

If we want to know these wave function at time t we have to apply the propagator associated 
to © to the initial state (fTT|) which we proceed to do in the next section. 



3 The wave function for the deuteron at t > 
in a suddenly apphed electrostatic field 

As we have obtained now the propagator Q, (|lflj) associated with the Lagrangian we can 
apply it to the state (fTT|) to get through equation (P) the state at a given time t > 0. 

The calculation involves integrals of exponentials of quadratic expressions in the variables f' 
and R' which can be carried out by completing squares in the exponents. We will calculate 
explicitly one example and then give the final result for ip{r,R,t) when t > 0. 

Remembering that for the x component of the relative vector r we have the propagator @, 
when we apply it to the part of the initial wave function exp(— ^wx^) of we have to evaluate 
the integral 

P AV3(--|-_y^%xp|i^ [(x2 + x'2) cota;t-2xx'(sin..t)-i]|e->"dx' 

= ^1/3 ( ^ \ ^ exp ( '—x^ cot ojt] r e-^^"+"^'dx', (12) 

V27risinu;t/ \ 2 J J _^ 
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where 



/9 = ^(1 — i cot wi), a = — ixu;(sinu;t) ^ 
and the real part of /3 is positive. From reference [5j we have that 
exp ( - + ax)dx = (vr//3)^ exp (a^/4/3) 



(13) 



(14) 



and from H13|) we obtain 



p = — (sma;t) exp(za;t), To ~ ^"^^ cotwt — -wx 

so that replacing in (|14|) and then in p2j) we obtain 
1 



a 



1 



^1/3 



exp 



-— 
2 



exp 



(15) 



which, as we should expect is, in our units, the phase term associated with the energy (i^/2) of 
the ground state of the oscillator. 

The same result holds for the y variable and both x, y can be replaced by x, y, as the 
electrostatic potential is applied only in the z direction. 

For the case when our initial wave function is A^l'^ exp(— u;z^/2), it is best to replace in it z 
by \z — (<f/\/2u;^)], as indicated in ©, and use for the variable z the propagator of the oscillator 
to which only a time dependent phase factor is multiplied as discussed after equation ©. 

With transition to the center of mass coordinates of the deuteron, denoted by capital let- 
ters, the propagator is just that of the free particle when we have X, Y so applying it to the 
corresponding part of the plane wave exp(ii('x)) exp(«i^yy) they will only give the phase factor 

exp(iK^X) exp ( - iK^t/2) , exp{iKyY) exp ( - iK^t/2) . 

For the Z center of mass coordinate we have to apply to exp(iKzZ) the propagator Q and the 
evaluation of the integral again follows procedures similar to those indicated in equations (|12j) 
to (HSJ. 

Combining then all our results we can say that for t > our wave function 'ip(r,R,t) will 
become 



i){r,R,t) = A[e a-""' e 



exp 



K R - -K'^t 



exp 



i£H \ 
'4uj^ J 



X exp 
X exp 



-ILUt 



^)^+2.3 
£t^ Z£t 



-iujt 



-(1 + COS^ U)t) 



+ 



i sin 2ijjt 



\- Kz 7= ^ 

12 "2V2 V2 



(16) 



4 The probability density for the deuteron 
in a suddenly applied electrostatic field 

If we turn to the "Diffraction in time" problem discussed in the Appendix, we note that to get 
information on its time dependent behavior we need not the wave function but its absolute value 
squared i.e. the probability density. Once we have this we can discuss the behavior in time with 
the help of the Cornu spiral [21. This of course is a general procedure in quantum mechanics 
and thus starting from the wave function in (fTB|l we need to write down \7p{r,R,t)\'^. 
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We see immediately that all part associated with the center of mass coordinates, i.e. capital 
X, Y, Z, disappears, as it would also disappear for the wave function at time t = given by 

For the relative coordinates we have to keep in the exponent of (|16() only the terms that are 
real multiplied by 2, so we get 

\'ilj(r,R,t)\'^ = ^^exp <^ -wr^ ^(1 - coswi) - 2coswt + (l + cos^wt)] } 

[ V2w 2uj-i^ ^-'J 



A"^ exp [— y^)] exp 

exp [ — uj[x'^ + y^)] exp 

A"^ exp [ — (jo(^x'^ + y^)] exp < 



-U! 



2£ S"^ 

z^H — ^— (1 — cosujt)z H j(l — coswi)^ 

V2u; 2a;'' 



-ijj 



-UJ 



Z H ^^(1 ~ COSLOt) 



n 2 



sin2(c^t/2) 



where used a trigonometric relation for the last expression. 

We have thus a very simple transient behavior that depends only on the relative coordinates 
and in which rr, y have the standard Gaussian behavior of range in our units, while z 

oscillates around this range with an amplitude (-v/2'?/'^^) sin^(u;t/2). The period of oscillation 
is 



(wr/2) = vr 



or 



T = (27r/a;) 



while the amplitude is (\/2<S/l<j^). 

Note that coordinates of proton and neutron, in the direction of the electrostatic field, are 
given respectively by 



so if we take our origin at the position of the center of mass, the distance between proton and 
neutron goes as \f2z and thus the deuteron is vibrating with an amplitude proportional to 
(2if /a;^) sin^(a;t/2) and it could radiate for strong electrostatic fields. 

Clearly for very strong electrostatic fields the deuteron could desintegrate and this information 
can be obtained also from the discussion of the classical version of the problem. 

The information available for propagators in reference [2] would allow us to discuss also the 
problem of a time dependent electrostatic field which reflect more the physical situation. The 
formulas we could use (6.2.42) is in p. 180 of reference [2]. 

For a constant electromagnetic field given by a vector potential A we could use formula 
(6.2.17) of p. 175 of reference 0. 

In general Feynman procedure gives a powerful tool for analyzing transient phenomena when 
external fields are applied to bound systems. 



5 General procedure for dealing with transient phenomena 
caused by a sudden perturbation 

While the propagator K(x,t;x') is very useful for the discussion of transient phenomena, it is 
not usually available for the general potential V+(a;). It is thus useful to develop an approximate 
procedure that is always available. 
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Denoting for H± the Hamiltonians associated with the potential V± we have that 

H_=T + V-, H+ = T + V+ = H_ + {V+-V-) 

where T is the kinetic energy. 

We shall assume that H- has only a discrete spectrum. In the case it also has a continuous 
spectrum we can reduce it to the previous case by inserting the whole system into a box. We 
shall denote the energies and eigenfunctions of H_ as 

H-4>n = En4>n, 

where 

Eq < El < E2 < ■ ■ ■ < En-l < En < En+l ■ ■ ■ ■ 

For Hj^ we are in search of a solution 'il){x, t) of the time dependent equation (in c.g.s. units) 

ih^ = = [H_ + {V+ - V.)]i; (17) 
with the initial value 

where v is some integer. 

We define the Laplace transform of ip{x,t) by 

/•oo 

■il}{x,s)= / e~^^ip{x,t)dt. 
Jo 

Applying the transform to equation ()17() we get 

ih J e'^^'^-^dt = ih J —{e~''^tp{x,t))dt + ihstp{x,s) 

= -ihipix, 0) + ihsi!{x, s) = [H^ + {V+ - V^)]ij{x, s). (18) 

Now the eigenstates (pn of H- are a complete orthonormal set so we can make the expansion 
00 

^(a;,s) = ^an{s)(pn{x) 

n=0 

SO that equation ((TH|) becomes 

-ihcpuix) +ihs^an{s)(l)n{x) = [H_ + {V+ - V-)]^an{s)(j)n{x). 

n n 

Multiplying by -0*, (a;) and integrating over the variables x we get 

-ih5n'u + ifisan'is) = En'an'{s) + '^{n\V+ - V-\n)an{s), (19) 

n 

where 

{n'\V+ - y_|n) = J cl)*^,{x)[V+{x) - V.{x)]Mx)dx. 
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The equation H19|) corresponds to an infinite set of linear algebraic equations for the coef- 
ficients a„(s). We can solve them with the usual assumption that after a given index N all 
coefficients vanish, and then we get the an{s) for n = 0, 1, . . . , from which we can write the 
tpix, s) as 



N 



n=0 



We want though to determine ip{x,t) for which we need the reciprocal of the Laplace trans- 
form 



ip{x,t) = 7^ / ij{x,s)e''^ds = — / eyi.^{-iEt/h)i){x,-iE/h)dE, 



AfC+ioo i-oo+ihc 
i}{x,s)e'*ds = 1^ 
c—ioo ZTTfl J_oo+ific 

where we replace the variable s by 

s = -i{E/h) 

with E being integrated along a line above all poles of the function '4){x, —iEh~^). Note that 
the contour can be completed with a circle from below and get essentially the residues at values 
of E given by the homogeneous set of linear equations 



{E - En')an' = ^(n'|V+ - F_|n)a„ 



or 



Ean' = y^^{n\H+\n)an 

n 

which essentially give us the energies associated with the levels of the Hamiltonian 

A Diffraction in time 

Units 

h = m = c = 1, 

u{l) = (h/mc), u{t) = {h/mc^), u{m) = m. 
Free particle propagator 

K{x,t;x',id) = (27rit)~^exp [{i/2t){x - xf]. 
The integral for diffraction in time is 

M{x,k,t) = / K{x,t; x ,0) ex.p{ikx')dx' , 



oo 



(27rit)-5 exp (^kx - ^fc^f^ | j exp + {kt - x)]"^^ dx' , 



where the last expression comes from completing the square for x' . 
We introduce the variable u by the definition 

(7r/2)5n = {2t)~^[x + {kt - x)] 
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from which it follows that for x' = 0, we have u = uq 
uo = {'7Tty^{kt - x), ^Jixjldu = {2t)~^dx 
and thus M{x,k,t) becomes 



exp 



i { kx — -k'^t 



i-uo 

exp(z7ru^/2)(in + / exp{iTru'^ /2)du 
■oo Jo 



But we have [3] 



° . 2 , ^ , (1 + ^) 

exp(i7ru /2)du = — - — . 



with the definition of Presnel integrals 



C{uq) = I cos{ttu /2)du, S{uo) 
Jo 

M{x,k,t) = e-'^'/^exp i\kx-\kH 



and 



exp(i7ru^/2)du = C{uq) + iS{uo) 



no 



sin(7rn^/2)du, 



\M{x,k,t)\' 



1 2 



+ 



n 2 



(20) 



for which the Cornu spiral can be used to show the difference between the quantum behavior 
given by (j2flj) and the classical one in which the probability density is 1 in the interval oo < x < vt 

and for vt < x < oo 
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